Experimental study of double-/? decay modes using a CdZnTe detector array 
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An array of sixteen 1 cm 3 CdZnTe semiconductor detectors was operated at the Gran Sasso 
Underground Laboratory (LNGS) to further investigate the feasibility of double-/? decay searches 
with such devices. As one of the double-/? decay experiments with the highest granularity the 4x4 
array accumulated an overall exposure of 18 kg days. The setup and performance of the array is 
described. Half-life limits for various double-/? decay modes of Cd, Zn and Te isotopes are obtained. 
No signal has been found, but several limits beyond 10 20 years have been performed. They are 
an order of magnitude better than those obtained with this technology before and comparable to 
most other experimental approaches for the isotopes under investigation. An improved limit for the 
/3+/EC decay of 120 Te is given. 



PACS numbers: 

I. INTRODUCTION 

Following the past 20 years vast progress has been 
made in unveiling the properties of the neutrinos. For 
decades neutrinos were thought to be massless, which no 
longer holds true: flavor oscillations found in the leptonic 
sector, studying neutrinos coming from the sun, the at- 
mosphere, high energy accelerators beams and nuclear 
power plants, are explained by a nonzero neutrino mass 
PIlElHISlHElll]. However no absolute mass scale 
can be fixed with experiments studying the oscillatory 
behaviour. To achieve this, one has to investigate on 
weak decays, such as beta decay or neutrinoless double- 
(3 decay (0z/2/3). The signature of the latter total lepton 
number violating process is the emission of two electrons 
with a sum energy corresponding to the Q-value of the 
nuclear transition. 

The 0z^2/?-decay is the gold plated process to distin- 
guish whether neutrinos are Majorana or Dirac particles. 
Furthermore, a match of helicities of the intermediate 
neutrino states is necessary, done in the easiest way by 
introducing a neutrino mass. This mass is linked with 
the experimental observable half-life via 

(T0; 2 ) = G^{Q,Z) \M% T - M°/\ 2 (^) 2 , (1) 

where (m ee ) is the effective Majorana neutrino mass, 
given by (m ee ) = \ ^ei m i| an d U e i as the correspond- 
ing matrix element in the leptonic PMNS mixing ma- 
trix, G 0,/ (Q, Z) is a phase space factor and M^ T - Mp 1 
describes the nuclear transition matrix element. Vari- 
ous stringent bounds on the half-life of several isotopes 
have been set. In addition, a potential evidence has 
been claimed in the 0^2/3-decay of 76 Ge with T°^ 2 = 



(2.23 ± 0.4) x 10 25 a at 90% C.L. (9]. For a recent re- 
view on double-/? decay see [10]. The COBRA experi- 
ment [TTJ uses CdZnTe (CZT) semiconductors to search 
for 0^2/3-decay. CZT contains nine double-/? emitters, 
five of which can decay via the emission of two electrons 
(in decreasing order of their Q-value these are 116 Cd 
, 130 Te , 70 Zn, 128 Te and 114 Cd), four of them can decay 
via double electron capture (EC) ( 108 Cd), a combination 
of EC and /?+ decay ( 64 Zn and 120 Te) or double /?+- 
decay ( 106 Cd). Evidently, every isotope which is able 
to decay with positron emission can also decay via dou- 
ble EC. The positron decays have a very nice potential 
signal of two (four) 511 keV annihilation gammas, but 
the phase space of the decay is reduced with respect to 
pure electron capture modes. There is a revived interest 
in these decay modes, as it has been shown that /3+/EC 
modes have an enhanced sensitivity to V + A interactions 
|12] and there might be a resonant enhancement in neu- 
trinoless double EC into an excited state of the daugh- 
ter nucleus if this state is degenerate with the ground 
state of the mother isotope [13] ■ The most promising 
isotope for the 0z^2/?-decay in COBRA is 116 Cd with a 
Q-value of 2809 keV resulting in a peak position beyond 
all gamma lines occurring from the natural decay chains 
of U and Th. There are a number of different experimen- 
tal approaches, also working on the aforesaid isotopes. 
The most advanced measurement for 130 Te is done by 
CUORICINO using cryogenic bolometers. The resulting 
half-life limit for the neutrinoless decay of 130 Te is given 
as T 1/2 > 3.0 x 10 24 yrs (90 % C.L.) [TJ]. Best limit 
for 116 Cd comes from enriched CdW04 scintillators |15| . 
These kind of crystals and also ZnW04 scintillators are 
used to set limits on the remaining Cd and Zn isotopes 
mentioned |16l IT?] . Recently a search for /3 + EC modes 
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and double EC modes of 120 Te has been performed using 
high purity germanium detectors. Results can be found 
in [18j. Most of the obtained limits are in the region of 
10 18 — 10 21 years. Results based on the semiconductor 
approach, namely the usage of CdZnTe detectors, were 
done with single detectors only [19] and on a small array 
of four detectors [20] . Here we report on results of a four 
times larger array consisting of a layer of 4 x 4 detectors. 



II. EXPERIMENTAL SETUP 

The data analyzed were taken with an array of 4 x 4 
CdZnTe crystals installed in the COBRA setup at Lab- 
oratori Nationali del Gran Sasso (LNGS). 

CoPlanar Grid (CPG)j21] type CdZnTe detectors were 
used. The detector crystals, manufactured by eV Prod- 
ucts [22] , were of cubic shape, measuring llxllxll mm 3 
with a total mass of 103.9 g. Detector anodes and cath- 
ode were gold-coated, the detector surfaces were passi- 
vated with a manufacturer-specific red lacquer to provide 
a highly resistive surface layer and to maintain long term 
stability. 

The array was embedded in a custom-made Delrin 
(POM) holder-structure, electrical contacting was done 
using thin Kapton- (Polyimide) Flex-PCB's and a low ac- 
tivity copper loaded glue which was developed in-house. 
The whole setup is surrounded by a 5 cm inner copper 
shield, followed by 20 cm of lead. A copper shield against 
electromagnetic interference (EMI) and a neutron ab- 
sorber consisting of 7 cm borated polyethylene completed 
the experimental setup. An overburden of 1400 m of rock 
at the LNGS location provided an excellent muon shield 
of 3500 m water equivalent. 

The electronic readout chain used was custom made 
specifically for the experiment. It consists of a 16 chan- 
nel preamplifier with integrated anode-grid subtraction 
circuit, low-noise supplies for grid-bias and high voltage 
and a VME-data acquisition system using custom built 
ADCs. 

We optimised all parameters relevant for detector per- 
formance — specifically cathode bias voltage, grid bias 
voltage and CPG subtraction circuit weighting factor — 
individually for each crystal. The whole setup was cal- 
ibrated regularly with 22 Na (511 and 1274.5 keV), 57 Co 
(122.1 keV) and 228 Th (238.6, 2614.5 keV). The 228 Th 
calibration is especially useful, since it provides several 
lines reaching from 84.4 keV up to the highest naturally 
occurring gamma energy of 2614.5 keV. 

During the measuring period of the experiment, most 
of the time 12 detectors were active, yielding a total ac- 
tive mass of 77.9 g. 

An in-depth description of the experimental setup, 
readout chain and detector calibration is given in |23j. 



Isotope and Decay 


Energy (MeV) 


FWHM (%) 


Effic. e(%) 


116 Cd to gs 


2.809 


3.5 - 8.8 


61.0 


130 Te to gs 


2.529 


3.6 - 8.9 


65.4 


130 Te to 536 keV 


1.993 


3.8 - 9.1 


58.2 


116 Cd to 1294 keV 


1.511 


4.1 - 9.4 


71.4 


130 Te to 1122 keV 


1.407 


4.2 - 9.5 


55.7 


116 Cd to 1757 keV 


1.048 


4.7 - 10.0 


60.5 


128 Te to gs 


0.868 


5.1 - 10.4 


92.1 


116 Cd to 2027 keV 


0.782 


5.3 - 10.6 


68.2 


130 Te to 1794keV 


0.735 


5.5 - 10.7 


63.0 


116 Cd to 2112 keV 


0.697 


5.6 - 10.9 


77.2 


116 Cd to 2225 keV 


0.584 


6.2 - 11.4 


77.1 


114 Cd to gs 


0.536 


6.5 - 11.7 


96.5 



TABLE I: 0v/3~ /3~ decay candidates and their decay energy 
together with the range of energy resolutions for all detec- 
tors (given as FWHM) at the decay energy and the average 
efficiency to detect the whole decay energy in one detector. 
Efficiencies vary typically within 1-2 % of the given value due 
to the different energy resolutions. 



III. DATA ACQUISITION 

The data analyzed in this paper can be divided into 
four periods. Small changes were applied during data 
taking which discriminate the periods from each other. 
After the first period, the whole setup was moved to a 
new location in the LNGS underground laboratory, but 
besides that, no changes were applied. For the third pe- 
riod, new preamplifiers were installed. During the fourth 
period, the detector channel assignment was changed. 
During periods one through three up to 12 crystals of the 
array were operational, while for the last period, only ten 
crystals were operational. 

Some cuts were applied to the data, to omit runs that 
have a high event rate due to e.g. microphonics result- 
ing in piezo-electric discharges. A run is equivalent to 
an hour of data taking, counted for each detector sepa- 
rately. The first cut applied did reject files written by 
the data acquisition software with a size of more then 
5 MB (a typical run has about 200 kB file size). In peri- 
ods one through three this applied to 2634 out of 43580 
runs. In the fourth period, 340 runs out of 25320 were re- 
jected. Furthermore runs were rejected on a detector by 
detector basis. Therefore, a Poissonian was fitted to the 
number of counts above 500 keV and runs with a count 
rate exceeding the Poissonian with a probability higher 
than 99% were rejected. Thus another 3157 runs out of 
65926 remaining runs were cut (4.8%). The spectrum of 

No further events 



all usable runs is shown in Figure 1 



were removed from the raw data. 

Variations observed in the resolutions could be at- 
tributed to the different preamplifiers used and the dif- 
ferent voltages applied in the different data taking peri- 
ods. It should be noted, that the detectors used during 



3 



>> 

(0 

> 



10 



c 

o 
o 



10" 1 



Summed spectrum, 18.0 kg days 



^'v'AViVj,.!,,,, 



500 1000 1500 2000 2500 3000 

Energy [keV] 



FIG. 1: The overall sum spectrum acquired with 12 detectors. 
Clearly visible is a low energy shoulder at 300 keV due to the 4- 
fold forbidden non-unique beta decay of Cd . Furthermore, 
two 7-lines at 352 keV 609 keV from the 238 U decay chain can 
be identified. 
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FIG. 3: Fit results for 116 Cd (2809 keV) and 130 Te (2529 keV) 
0v/3~ f3~ at Figure a) (combined fit for both isotopes), and 
106 Cd EC/EC decays at Figure b). For the /?"/?" isotopes, 
a Gaussian signal was fitted together with a background 
parametrization. For the /3 + /3 + modes, simulated data sets 
were used to describe the signal. 




FIG. 2: Simulated energy spectrum of 120 Te Qvf3 + fEC to 
ground state as an example of the complexity to be expected 
in such decay modes. The decays were simulated for each 
detector individually assuming ideal resolution, but taking 
into account the setup geometry of the whole detector layer. 
Afterwards, the spectra were convoluted with the resolution 
of the given detector for each run individually. The spectrum 
is described in the text. 



these measurements do not have the best energy reso- 
lution possible for CZT CPG semiconductor detectors. 
A better energy resolution would be desirable, but not 
essential at the moment, since the background was not 
known and therefore cheaper detectors were used. 



IV. DATA ANALYSIS 

From the cleaned data sample of 18 kg days, upper lim- 
its on the number of signal events were extracted using 
a maximum likelihood fit. Therefore, the data was his- 
togrammed for each detector. For the fit, each bin of 
a histogram was modeled as a Poisson random variable 
depending on the signal and background rate. The com- 
plete likelihood function, deriving from the Poissonian 
probability distribution, for all detectors and all runs dur- 
ing the measurements is given by 

ln(L) = S 2 ^ Bin ' ln Fm ^ ~ Fb ^ ^ 

Det. Run Bin 

with the number -Dein of observed events in the given 
histogram bin and the fitted number Fein of counts de- 
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termined from the fit function. This approach takes 
into account the energy resolution and background level 
of each detector individually by the parameters of i<Bm 
and also allows the consideration of changes of resolution 
and background levels during different runs. The log- 
likelihood function was minimized with the MINUIT |24| 
package. 

The fit model for Fbui consists of a background 
parametrization together with a signal parametriza- 
tion. No background was subtracted from the ana- 
lyzed data sample. The background model contains two 
exponentials, one describing the higher energetic spec- 
trum above 700 keV, the other describing the spectrum 
above 400 keV. Additionally, a Gaussian is fitted to the 
609.3 keV 214 Bi line. Below 400 keV, the spectrum is 
mainly dominated by the fourfold forbidden /3-decay of 
113 Cd |23 and low-energetic 7 lines like the 351.9 keV 
214 Pb line. As the lowest peak investigated is at 534 keV, 
data below 400 keV were not taken into account for the 
analysis. 

For the signal model, different approaches for 0^/3~/3~ 
and Qv(3 + (3 + decays were taken. The summed energy of 
the two electrons of a Qvf3~ (3~ decay results in a Gaus- 
sian signal peak with a mean at the decay energy E pea k- 
Such a Gaussian was taken as signal model for the (3~ (3~ 
decays. To ensure a proper fit of the background, a fit 
range of at least three peak widths (FWHM) around the 
expected peak energy was chosen. Due to the close prox- 
imity of the expected signal peaks of 116 Cd and 130 Te as 
well as the 70 Zn and 128 Te ground state transitions, a 
combined fit of both peaks was made for these signals. 
To obtain the efficiency e for observing the full decay 
energy, Monte Carlo simulations taking into account the 
escape probabilities of the electrons and also of result- 
ing gammas were done for each detector. The expected 
Qvj3 + f3 + , namely double positron transition /3 + /3 + , sin- 
gle electron capture /3+/EC and double electron capture 
EC/EC spectra are more complex than the f3~ (3~ ones. 
They normally contain several lines deriving e.g. from 
the escape of 511 keV annihilation gammas from the de- 
tectors, bu t do not h ave a dominating peak. An example 
is given in |Figure~2 with the simulation of 0^/3 + /EC de- 
cays of 120 Te. There, besides a peak at the full Q-value of 
1722 keV and the mentioned escape peaks, also satellite 
lines 22 keV below these lines deriving from additional es- 
cape of a K a X-ray and a peak at 511 keV deriving from 
gammas produced in annihilations of positrons from the 
/3 + decays in other detectors of the 16er layer are visible. 

To cope with these complex signal distributions, an 
alternate approach for the search for these signals was 
taken. The spectra of the 0vf] + f3 + decay modes were 
simulated for each detector taking into account the mass 
of each detector and its position inside the layer. The 
resulting spectrum was convolved with the energy res- 
olution of the given detector, normalised to unity and 
taken as signal model. The fit range for every isotope 
includes the main lines of the given decay. 

Results of the fit of the background and the signal 



model for the 0^/3~/3~ decays of 116 Cd an d 130 Te a nd 
for the (3 + /EC decay of 106 Cd are shown in Figure 3 

The simulation of /3 + /3 + and f}~ (3~ signals in the detec- 
tors of the experimental setup were done with a GEANT4 
based Monte Carlo simulation. The Fortran DecayO [25 
code was used as event generator. 

A 90% confidence level upper limit on the fitted counts 
was derived from the shape of the log-likelihood func- 
tion. It is defined as the point where the (negative) log- 
likelihood function rises by a value of 1,28 above its min- 
imum. The shape of the log-like lihood fun ction for the 

In case the 



106 



Cd Oj/EC/EC fit is shown in Figure 4 



count rate obtained by the fit is positive, the 90% C.L. 
upper limit 7V s i g on signal events is calculated by sum- 
ming the fit result and the associated 90% C.L. error. 
In case of it being negative, a conservative approach by 
only using the 90% C.L. error has been applied. With 
the obtained upper limit, a 90% lower half-life limit can 
be calculated via 



Tl/2 > 



ln(2) • AT i6 



(3) 



where Aj SO is the number of source atoms under study, 
tiivo is the lifetime of the experiment and e is the above 
mentioned simulated efficiency o f observ ing a decay in 
the experimental setup as given in Table I For the f3 + f3 + 
decays, the whole simulated spectrum was used for the 
extraction of excluded events, thus e = 1 holds for these 
calculations. 

The Cd and Zn content of the Cdo.9Zno.1Te is only 
known to lie in a range of 7% to 11% for Zn and ac- 
cordingly 89% to 93% for Cd. Therefore, conservatively 
a zinc content of 7% was assumed when calculating the 
limits for Zn isotopes and 11% zinc content was taken for 
the Cd isotopes. 



RESULTS 



The obtained results for j3_ 
blellland for /3+/3+ 



/? decays are listed in Ta- 



Table III Due to the low source 



mass of 77.9 g not all results can compete with the results 
of world leading large scale experiments. However, for- 
mer COBRA limits could be significantly improved, six 
of them by more than an order of magnitude and several 
have passed the 10 20 yr boundary. In addition, various 
decay channels are within an order of magnitude of the 
world leading limits. 

Due to significant improvements possible in the near 
future there is a good chance to surpass some of the ex- 
isting limits. The installation of a nitrogen atmosphere 
including a radon trap combined with the replacement of 
the red passivation lacquer indicates a reduction of back- 
ground by an order of magnitude. The installation of 48 
more detectors will add another factor three in source 
strength and also significantly increase the sensitivity for 
several decay modes especially those involving photons. 
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FIG. 4: The shape of the negative log-likelihood function for 
the EC/EC fit of 106 Cd. The 90% upper confidence level is 
determined by the point where — ln(L) rises 1,28 over the 
minimal value. 



Isotope and Decay 


Fit Range 
(MeV) 


T1/2 limit (years) 


This work 


Measured 


116 Cd to gs 
130 Te to gs 


2.2-3.2 
2.2-3.2 


9.4 x 10 19 
5.0 x 10 20 


1.7 x 10 23 fj5] 
3.0 x 10 24 ES] 


130 Te to 536 keV 


1.7-2.3 


3.5 x 10 20 


9.7 x 10 22 E3 


116 Cd to 1294 keV 


1.2-1.8 


5.0 x 10 19 


2.9 x 10 22 fTS] 


116 Cd to 1757 keV 


0.9-1.3 


4.2 x 10 19 


1.4 x 10 22 fTS] 


128 Te to gs 


0.6-1.3 


1.7 x 10 20 


1.1 x 10 23 [H] 


116 Cd to 2027 keV 


0.5-1.2 


2.8 x 10 19 


2.1 x 10 21 fJS] 


116 Cd to 2112 keV 


0.5-1.0 


4.7 x 10 19 


6.0 x 10 21 rrs] 


116 Cd to 2225 keV 


0.5-1.0 


2.1 x 10 19 


1.0 x 10 20 t (30] 


130 Te to 1794keV 


0.5-1.2 


1.9 x 10 20 


2.3 x 10 21 ED 


130 Te to 1122 keV 


1.1-1.7 


1.2 x 10 20 


2.7 x 10 21 OLD 


114 Cd to gs 


0.4-1.0 


2.0 x 10 20 


1.1 x 10 21 fj7] 



TABLE II: Results for fits of Qvf3~ (3~ decays on 90% C.L. and 
fit ranges. Entries marked with ' are on 68% C.L. Isotopes 
not separated by a horizontal line were fitted together. 



Isotope and Decay 


Fit Range 
(MeV) 


T1/2 limit (years) 


This work 


Measured 


64 Zn /3+EC to gs 


0.5-1.1 


1.1 x 10 1S 


4.3 x 10 20 fj6] 


64 Zn 2EC to gs 


0.5-1.3 


3.3 x 10 17 


1.1 x 10 20 [E] 


120 Te (3+EC to gs 


1.0-2.0 


4.1 x 10 17 


1.9 x io 17 in] 


120 Te 2EC 


0.8-2.0 


2.4 x 10 16 


6.0 x 10 17 [18] 


120 Te 2EC to 1171 keV 


0.6-2.0 


1.8 x 10 16 


6.0 x 10 17 [18] 


106 Cd (3+ (3+ to gs. 


0.5-2.0 


2.7 x 10 18 


2.4 x 10 20 (32] 


106 Cd (3+EC to gs 


1.5-3.0 


4.7 x 10 18 


3.7 x 10 20 (32] 


106 Cd 2 EC to gs 


2.0-3.0 


1.6 x 10 17 


3.5 x 10 18 (33] 


106 Cd f3 + /3+ to 512 keV 


0.6-1.5 


9.4 x 10 17 


1.6 x 10 20 (32] 


106 Cd /3+EC to 512 keV 


0.8-2.0 


4.6 x 10 18 


2.6 x 10 20 (32] 



TABLE III: Results obtained for Quf3 + (3 + modes on 90% C.L. 
and fit range. 



detectors. For such experiments a low background rate 
in the peak region and a good energy resolution are cru- 
cial ingredients. To further improve this approach, CZT 
semiconductor detectors have been installed and oper- 
ated in a low background environment deep underground. 

An 18 kg day data set was collected with up to 12 
1 cm 3 crystals running simultaneously. It was analyzed 
to determine half life limits on a number of neutrino- 
less double-/3 decay modes for seven different isotopes. 
Nearly all of the derived limits improved by roughly an 
order of magnitude compared to previous results |20j. 
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Furthermore, a potential coincidence analysis might al- 
low for even better background reduction. 



VI. SUMMARY 



The newly planned double-/3 decay experiment CO- 
BRA consists of a large amount of CZT semiconductor 
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